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Particle motion in a cylindrical multiple-cusp magnetic 
field configuration is shown to be highly (though not com- 
pletely) chaotic, as expected by analogy with the Sinai bil- 
liard. This provides a collisionless, linear mechanism for phase 
randomization during monochromatic wave heating. A gen- 
eral quasilinear theory of collisionless energy diffusion is de- 
veloped for particles with a Hamiltonian of the form Hq + Hi, 
motion in the unperturbed Hamiltonian Ho being assumed 
chaotic, while the perturbation Hi can be coherent (i.e. not 
stochastic). For the multicusp geometry, two heating mecha- 
nisms are identified — cyclotron resonance heating of parti- 
cles temporarily mirror-trapped in the cusps, and nonresonant 
heating of nonadiabatically reflected particles (the majority). 
An analytically solvable model leads to an expression for a 
transit-time correction factor, exponentially decreasing with 
increasing frequency. The theory is illustrated using the ge- 
ometry of a typical laboratory experiment. 
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I. INTRODUCTION 

The quasilinear diffusion equation in its original form 
(QLT1) was a Fokker-Planck equation describing the 
velocity-space diffusion of particles due to random scat- 
tering by waves. In the absence of the waves the plasma 
was assumed to be infinite and homogeneous so that the 
unperturbed motion was rectilinear. The diffusion equa- 
tion was derived from the Vlasov equation by solv- 
ing for the perturbed part of the distribution function, in 
the linear approximation and assuming the unperturbed 
distribution function to be essentially constant over the 
timescale of wave-particle interaction, then substituting 
back into the full Vlasov equation and averaging over 
position (or, equivalently, over the random phases of the 
waves). 

To put this formalism in a more general perspective, it 
is advantageous to cast it in Hamiltonian form, with the 
particle motion in the absence of waves being described 
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by an unperturbed Hamiltonian Hq. The total single- 
particle Hamiltonian is then Hq + ei?i, with the waves 
being incorporated in the perturbation Hamiltonian H±. 
In QLT1, all stochasticity comes from the assumption of 
random phases in the assumed broad spectrum of waves 
in Hi . The smallness parameter e expresses the assump- 
tion that the amplitudes of individual waves are small, 
so that the short-time-scale perturbations to the distri- 
bution function can be derived using linear, O(e), theory 
(hence the terminology "quasilinear"). The only non- 
linear effect is the long-time diffusion described by the 
diffusion coefficient, which is 0(e 2 ) . 

It is not really necessary to assume an infinite, homo- 
geneous plasma. The essence of QLT1 is the assumption 
that Hq is integrable, so that a canonical transformation 
to action-angle coordinates exists. Then the unperturbed 
motion in angle space is rectilinear, just as that in coor- 
dinate space in the case of a homogeneous medium. This 
action-angle generalization was carried out by Kaufman 
B and by Hazeltine et al. || to derive a quasilinear dif- 
fusion equation (in action space) for wave-particle scat- 
tering in axisymmetric toroidal magnetic confinement ge- 
ometries (e.g. tokamaks). The formalism has been used, 
for example ||, to investigate the effect of a sheared ra- 
dial electric field on anomalous transport in a tokamak. 

With the development of the theory of Hamiltonian 
chaos it has come to be realized that a quasilinear dif- 
fusion equation can also be derived in cases where Hi 
represents the effect of a coherent wave, provided the in- 
teraction of Hq and Hi produces a chaotic motion. We 
call this form of quasilinear theory QLT2, and it is very 
useful in the theory of radio- frequency (rf ) and microwave 
heating of plasmas because this is typically done with 
coherent waves. 

Again, Hq is assumed integrable so that an action- 
angle transformation exists. In these theories the pertur- 
bation is still the source of chaos ("intrinsic stochastic- 
ity"), which causes the action variables, constructed in 
the integration of Hq, to perform the random walk de- 
scribed by the quasilinear diffusion equation. Assuming 
the Hessian matrix d 2 Ho/dpidpj to be nonsingular, the 
coherent perturbation must exceed a certain amplitude 
for global resonance overlap, and hence chaos, to occur 
P-|lO|. Thus, paradoxically, a criterion for QLT2 to ap- 
ply is that the system be sufficiently nonlinear! 
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In the present paper we examine a third form of quasi- 
linear theory, which we call QLT3. This case is the com- 
plete obverse of the original quasilinear diffusion prob- 
lem: we now assume the unperturbed Hamiltonian Hq 
to be completely nonintegrable, giving rise to strongly 
chaotic motion in the unperturbed system. The chaotic 
dynamics of the unperturbed Hamiltonian system then 
provides randomization and allows the application of a 
quasilinear formalism to derive the diffusion equation for 
the distribution function. 

Of course, since Hq is not integrable, an action-angle 
representation does not exist. In fact we assume Hq to be 
so strongly chaotic that the unperturbed motion covers 
essentially the entire energy surface Hq = £ ergodically, 
except as restricted by integrals of the motion associated 
with any continuous symmetries of Hq. The goal of this 
paper is to determine the diffusion in £ caused by the 
time-dependent perturbation H\. 

Since the unperturbed motion provides the source of 
stochasticity (with no threshold), we can, as in QLT2, 
assume the perturbation to be coherent. Thus the the- 
ory is applicable to wave heating of plasmas in strongly 
nonintegrable magnetic confinement geometries. 

The assumption of uncorrelated gyrophase in succes- 
sive passes through the resonance region is often made in 
the derivation of a quasilinear diffusion equation to de- 
scribe cyclotron resonance heating of magnetically con- 
fined plasmas. However, in simple confinement geome- 
tries, such as magnetic mirrors, Hq is essentially inte- 
grable owing to the existence of the adiabatic invariant fi 
and another integral due to symmetry, or a second adi- 
abatic invariant. Lichtenberg and Lieberman |llj have 
analyzed collisionlcss heating in such systems using area- 
preserving maps, and find the random phase assump- 
tion to be valid only well beyond the nonlinear threshold 
where the last invariant circle is destroyed and chaotic 
motion becomes global. In our nomenclature, quasilin- 
ear diffusion in these systems is an example of QLT2, not 
QLT3. 

On the other hand, in systems with a null in the mag- 
netic field /i is not globally conserved and the situation 
is rather different from that in the much-studied mirror 
systems. For instance, Yoshida et al. |[2| have recently 
studied rf heating in a two-dimensional slab model with 
a neutral line, where fi conservation is broken. They find 
heating due to the onset of chaos at finite amplitude of 
the perturbing field. However, their unperturbed system 
has two symmetry directions, and thus their Hq is inte- 
grable, despite the breaking of the adiabatic invariant. 
Thus their model must be classified as a QLT2 case also. 

In systems with a null in the magnetic field and only 
one symmetry direction, however, Hq is not in general 
integrable. An important class of such systems are the 
magnetic cusp confinement geometries, which are much 
used in low-temperature plasma physics due to the ease 
with which they can be created with arrays of perma- 
nent magnets |l3|] jl4|, pp. 146-150]. In this paper we 
give evidence that they fulfill the criterion for systems of 



type QLT3 in that their unperturbed dynamics is almost 
completely chaotic. 

We know from the work of Sinai |l5|,[l6| that particle 
motion on a billiard table with a convex boundary is a 
strongly chaotic system due to the defocussing effect of 
each collision with this boundary. In fact Sinai showed 
the motion to be strongly mixing, so that ergodic theory 
could be applied. This suggests that cusp confinement 
systems, where magnetic fields lie on surfaces that are 
convex toward the plasma, are strongly chaotic. Indeed 
the four-cusp Hamiltonian H = \{p\ + x 2 y 2 ) was at 
one time conjectured to be completely chaotic like a Sinai 
billiard. Although Dahlqvist and Russberg [jl7| disproved 
this conjecture by finding a stable periodic orbit, they 
found that the island of stability surrounding this orbit 
was extremely small, so for practical purposes one can 
assume that the energetically accessible phase space is 
covered ergodically in this system. 

In this paper we study a magnetic field configura- 
tion consisting of a "picket fence" of infinite linear mag- 
netic dipoles, producing multiple line cusps. This can 
be regarded as a simplified model for a low-temperature 
plasma source, created using arrays of permanent mag- 
nets [ff8 19 1 . It is also a rather simple model, in which 
the unperturbed dynamics can be simplified analytically 
to a considerable extent by the use of complex variable 
theory. 

In the limit of a large number of dipoles the interior 
of the plasma is essentially free of magnetic field and the 
unperturbed motion is rectilinear in the interior region, 
while near the edge of the confinement region, a par- 
ticle can be reflected over a range of angles. Thus we 
might expect the configuration to approximate a chaotic 
billiard problem. (In the original Sinai problem the con- 
vex boundary was in the interior of the domain, whereas 
the billiard analog of the present example has an outer 
boundary that is convex except for cusps, like the "bow- 
tie" billiard shown in Fig. 7.24(e) of Ref. @.) 

The problem is also related to a model originally pro- 
posed by Fermi pc|-|22]| for explaining the acceleration of 
cosmic rays to the extraordinarily high energies observed. 
In the Fermi model the cosmic ray particles move recti- 
lincarly except during occasional collisions with moving 
magnetized clouds of gas, which cause diffusion in energy 
space. The present model includes both the possibility 
of cyclotron-resonance heating in the mirror-like cusp re- 
gions and nonresonant heating of particles reflected with- 
out penetrating deeply into the cusps. The latter case is 
much closer to the original Fermi problem and is the main 
focus of the paper. 

In Sec. H we introduce the confinement geo met ry and 
unperturbed Hamiltonian in detail, and in Sec. Ill we an- 
alyze the dynamics of this system and show it is indeed 
strongly chaotic for the class of particles ( "free particles" ) 
that traver se the central low- field region. However, in 
Sec. HID 3 we produce a counter example to any conjec- 
ture that the motion is totally chaotic by finding a stable 
periodic orbit. 
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In Sec. IV we introduce the wave-particle interaction. 
The general quasilinear diffusion equation is derived in 
Sec. M. An analytically solvable one-dimensional model 



of the magnetic field is used in Sec. VI to estimate heat- 
ing of nonadiabatically reflected particles. The result is 
of the form expected from simple Fermi acceleration the- 
ory multiplied by a transit-time reduction factor that be- 
comes exponentially small when the transit time is much 
longer than the period of the applied field. The the- 
ory is discussed using typical parameters for permanent- 
magnet confinement experiments in Sec. VII. 



II. UNPERTURBED HAMILTONIAN 
A. Two-dimensional magnetic Hamiltonian 

The behavior of a sufficiently dilute plasma can be an- 
alyzed on the basis of single-particle motion in magnetic 
and electric fields made up of an externally imposed com- 
ponent and an internally generated component produced 
by the collective currents and charges from the combined 
effect of many otherwise noninteracting particles. In this 
paper we suppose that the self-consistent component is 
negligible and analyze single-particle motion in an im- 
posed magnetic field. 

Consider the motion of a particle of charge q and mass 
m in a straight, infinitely long magnetic confinement sys- 
tem with vector potential A = ip(x, y)e z , where x, y and 
z are Cartesian coordinates and e z is the unit vector 
in the z-direction. The magnetic field is B x = dtp/dy, 
By = —dip/dx, B z — 0, so contours of the flux function 
ip(x, y) in any plane z = const define magnetic field lines. 

The Hamiltonian is (see e.g. Q ) 



H a = 



P 2 

fx 

2m 2m 



P 2 v , \Pz-q^{x,y)f 



2m 



(1) 



where pi (i € {x,y,z}) are the canonical momenta, 
Hamilton's equations of motion being Xi — dH^/dpi, 
pi = -dH /dx l . 

B. Multicusp flux function 

Assuming no currents present in the plasma, -0 obeys 
Laplace's equation. It is a standard result that two- 
dimensional solutions of Laplace's equation can be con- 
structed by taking the real or imaginary part of any an- 
alytic function of £ = x + iy p4| . Thus we write 



Mx,y) =Re*(C) , 



(2) 



where ty(x,y) is the complex flux function. 

For instance, a line current (line magnetic monopole) 
at Co is represented by the real (imaginary) part of ln(£ — 
Co). 



Although a magnetic monopole cannot be realized in 
nature, a magnet can be modeled as a superposition of 
magnetic dipoles. We consider a magnet that is long 
in the z-direction, thin in the x- and y-directions, and 
magnetized in the ^-direction so that it can be modeled 
by a line magnetic dipole. Such a linear dipole can be 
produced by differentiation of Im ln(£ — (o) with respect 
to Xq. 

Superimposing the flux functions for 2n linear magnets 
of strength alternately +K and —K lying in a circular 
cylinder of radius a about the z-axis we find ifi for a 
circular multidipole magnetic confinement system, 



ip = Kim 



n-1 

E 



An 



An' -I 



2nK 



-Re 



An 1 

au n 



au r , 



(3) 



Here u n = exp(7ri/2n) is the 4nth root of unity. The 
equivalence of the first and second forms can be verified 
by showing that they have the same poles and residues 
and the same value at Q — 0. 

Thus, comparing with Eq. (||), we see that the complex 
stream function for a circular dipole picket fence is given 
by 



2nK 



a 

nK 
a 



sech 



In 



(4) 



It is clear from the first form that ^ has poles at the 
2nth roots of — a 2n . In terms of polar coordinates r and 
such that £ = rexpi#, the poles are at r = a, 9 — 
n/2n + 2/7r, I = 0, 1, . . . , 2n — 1. Contours of (lines 
of force) are shown in Fig. ^ for the case n = 6, with 
distances measured in units of a/n. 




3 4 
X 



FIG. 1. Contours of the flux function for a typical cylindri- 
cal multicusp geometry produced by 12 line dipoles of alter- 
nating polarity. 



C. Near field and nondimensionalization 

We now consider the behavior of in the region be- 
tween two magnets, which for definitcness we take to be 



3 



those immediately above and below the positive real x- 
axis. We expand about the point £ = a, on the circle 
on which the magnets are located, by setting £ = a + £, 
where ^ = x — a -\- iy. Assuming |£| <C a, we see from 
Eq. (g) that 



sech I — 

a \ a 



(5) 



Thus, on the a;-axis ip peaks at x = a and it decays 
rapidly on either side. Furthermore, we see that the 
scale length of the x and y variation is not a but a/n. 
In Sec. Ill D 1 we shall encounter Eq. (J^) again in the 
context of the asymptotic limit n — > oo. 

Expanding W to second order in n£/a we find ip = 
Re\E' w f/'xil — ^ 2 [(a;— a) 2 — y 2 ]/2a 2 } in the neighborhood 
of (a, 0), where 



?/>x = nK/a 



(6) 



is the value of t/> at this saddle point, the location of a 
magnetic-field null. It is also useful to define a typical 
magnetic field Bq in the strong-field region via the rela- 
tion Bq = nipx/a- 

We term the energy of a particle with momentum p = 
located at this saddle point the escape energy 



f = 



<?V 2 



|_ 1 fqnKV 
2m 2m \ a J 



(7) 



In all numerical work and figures in this paper we 
nondimensionalize by measuring distance in units of a/n, 
mass in units of m, time in units of the typical inverse 
angular cyclotron frequency 



\q\B n 2 Kq ' 



(8) 



and charge in units of \q\. In these units a — n, m 
\q\ = r x = ipx = B Q = 1, and £ csc = §. 



B. Integrability and effective potential 

Because Hq is independent of z, p z is a constant of 
the motion. Also, Ho itself is a constant of the motion, 
Hq = £, where the constant £ is the total energy. How- 
ever, the absence of a third integral of the motion means 
the system is not integrable. Thus we must resort to nu- 
merical integration to investigate the nature of the un- 
perturbed orbits. 

Before proceeding to a discussion of the numerical re- 
sults, however, we observe that some qualitative under- 
standing of the motion can be found by determining 
the bounds of the motion implied by the constancy of 
H . Inspecting Eq. (|j) we see that the term V e s = 
(p z — qip) 2 /2m in H acts as an effective potential in which 
the particles move. The motion in the (x, y)-plane is thus 
bounded by the curves V e g (x, y) = £ . Note that V e g > 0, 
with equality occurring on the contours ip = p z /q- Also, 
since ip vanishes at the origin, V e g = p\/2m there. 

In the case p z — 0, the curves V e g(x,y) = const are 
just level curves of \tp\. There are thus two cases: 

1. Perfect confinement 

£ < £ csc , the curve \ip\ = (2m£) 1 ^ 2 /\q\ encloses the 
origin (though it has cusps at the magnets) and 
there is thus no leakage of particles with p z — 
through the dipole picket fence. 

2. Partial confinement 

For £ > £ esc , the curves = (2m£) 1 / 2 /\q\ are 
disjoint and thus particles can escape through the 
"mountain passes" (see Fig. |l|) between the mag- 
nets. 

In either case the particles are free to traverse a large 
region including the origin, like particles rolling on a bil- 
liard table [in case (ii) it is a billiard table with pockets] . 
We refer to particles on such orbits as free particles, to 
be contrasted with the trapped particles discussed in the 
next section. 



III. UNPERTURBED PARTICLE DYNAMICS 

A. Dynamics in complex notation 

For a given p z , the dynamics is a two-degree-of- freedom 
autonomous Hamiltonian system which can be compactly 
written using complex notation as 



c 



1 



■pc 



(9) 



PC = — [Pz 
m 



gRett(C)] [*'(C)r, 



where the prime on 4" means the derivative with respect 
to its argument, * means complex conjugate, and pq = 
Px + iPy 





FIG. 2. A typical trapped-particle orbit in the magnetic 
field shown in Fig. The thick lines show the boundaries of 
the energetically accessible regions. 
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C. Trapped-particles 

For p z ^ a new class of orbit arises, the trapped par- 
ticles. This occurs when the energy is less than the value 
of V e s at the origin, i.e. when £ < p\ /2m, for then the 
low-magnetic- field c entral region is energetically inacces- 
sible (see Sec. [IIP ) and the particle must be confined in 
the edge region near the magnets, as illustrated in Fig. ||. 

Because deeply trapped particles are always in a re- 
gion of strong ma gne tic field, their dynamics can be an- 
alyzed (see, e.g., [£3|, pp. 21-22]) by decomposing the 
motion into that of the guiding center, with velocity 

v = W|jB/_B + drifts, and a gyro motion with velocity 

vi in the plane locally orthogonal to B. The adiabatic 
invariant 



2B 



(10) 



is conserved to high degree of approximation, providing 
an approximate third integral of the motion. The un- 
perturbed dynamics of this class of orbit is thus quasi- 
integrable, not chaotic, implying that we must use quasi- 
linear theory of type QLT2 to derive a diffusion equation 
(i.e. heating will occur only beyond a nonlinear ampli- 
tude threshold). Cyclotron resonance heating in mirror 
geometries has been much discussed in the literature fl4| , 
pp. 413-422] and we shall not discuss the heating of the 
trapped particles in this paper. 



D. Free particles 

For a given energy £ and conserved momentum 
p z , the energetically accessible region is the set of 
points (x, y) for which there exist p x and p y such that 
H (x,y,p x ,py,p z ) — £. From Eq. (|l|) we get the inequal- 
ity 



[p z - qip(x,y)] 2 < 2m£ 



(11) 



We define the free particles as those for which the ori- 
gin is energetically accessible. Thus, since ^(0,0) = 0, 
Eq. (O) implies that free particles are those for which 



P 



I < 2m£ 



(12) 



The total region accessible to free particles is the set 
of (x, y) for which the ranges of p z defined by Eq. ( p"l] ) 
and Eq. ( [12] ) are not disjoint. This gives the condition 



\q>l>(x,v)\ <2(2m£) 1 ' 2 



(13) 



This inequality being satisfied, the intersection of the 
ranges defined by Eq. and Eq. (|l|) is [qip - 
(2ms) 1 / 2 , (2m£) 1 / 2 ] for qijj > 0, or [-{2m£) 1 / 2 , + 
{2m£) 1 / 2 ] for # < 0. 

A typical trajectory for a free particle with energy well 
below the escape energy is shown in Fig. 0. In this case 



p z = 0, and the initial conditions are xq = yo = 0, p x o 
0.04, p yQ = 0.05, giving £ = 0.00205. 



I\4 
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FIG. 3. A typical free-particle orbit in the magnetic field 
shown in Fig. m, with initial conditions as described in the 
text. 

We see that the energetically accessible region is 
bounded by a series of curves (shown as thick lines) join- 
ing in cusps at the magnets. These bounding curves 
are convex toward the confinement region, except at the 
cusps, where the magnetic field goes to infinity so that no 
particle can penetrate. Also we see that the motion of the 
particle well away from the bounding curves is approxi- 
mately rectilinear so that the system does indeed appear 
like a physical realization of a Sinai billiard system |l5| ] . 

Although the orbit in Fig. || looks chaotic, a better 
test for chaos is to do a Poincare surface of section punc- 
ture plot, as shown in Fig. ^ for a particle with p z = 
and energy £ = 0.0017 started near the periodic orbit 
described in Sec. [II D 3| . 
tion is x > 0, y 
operation £ 



The Poincare surface of sec- 
and its images under the symmetry 
exp(i7r/6)C- Dots indicate both upward- 
and downward-going passes of the orbit. It is seen that 
the orbit appears to fill the energetically accessible phase 
space ergodically, indicating strong chaos. 



Px 
0.06 »_ 



0.04 
0.02 

-0.02 
-0.04 
-0.06 f 



FIG. 4. Intersections of an orbit with the Poincare surface 
of section y — described in the text. 

However, the interaction with the field near the bound- 
ing curve is not specular reflection with a zero radius of 
curvature, so our magnetic cusp confinement system is 
not precisely analogous to that studied by Sinai. In fact, 
studying Fig. |§| we see that there are two qualitatively 
different kinds of reflection event — approximately spec- 
ular reflection with a small-but-finite radius of curvature, 
and "sticky reflections" near the cusps, where the parti- 
cle is temporarily trapped in a one-sided mirror field and 
oscillates several times before reflection. As explained be- 
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low, we call these nonadiabatic and adiabatic reflections, 
respectively. 

1. Scattering analysis of reflection, large n 

In order to make the reflection process a precisely de- 
fined event we go to the large- n limit, in which the spac- 
ing between the magnets, ira/n, and the scale length of 
magnetic field variation, a/n, become small compared 
with the radius a. Since we are interested in dynam- 
ics near the wall of magnets, we shift the origin to the 
sad dle p oint at £ = a by defining £ = £ — a, just as in 
Sec. EIC. Again Eq. (J|) applies at leading order, but this 
time its region of validity extends beyond the region of 
the magnetic null to include the high-field regions near 
the magnets (the only requirement being |£| -C a). 

We now take Eq. (|^) to be the exact complex flux func- 
tion for the model problem of an infinite line of magnets 
(treating n as an arbitrary parameter, which is scale d 
out in the nondimensionalization defined in Sec. Ill B| ) . 
A reflection event is now precisely defined as a scattering 
process, in which a particle impinges from Re£ = — oo 
(where the orbit is asymptotically a straight line), reflects 
off the magnetic field of the magnets, then retreats back 
toward Re£ = — oo. 




FIG. 5. A set of nine orbits, incident in the direction nor- 
mal to the line of magnets, with different "impact parameter" 
Ho. The thick line shows the boundary of the energetically ac- 
cessible region. The origin of the £-axis has been shifted to 
lie on the line of magnets. 

In Fig. H we illustrate some typical reflection events 
for particles with initial momentum pq x = 0.06, p y o = 
0, p z — 0, giving energy £ = 0.0018. For small initial 
values of y, yo, the reflection is approximately specular, 
but as yo approaches ir/2 (the height of the first magnet 
above the a:- axis), the orbit undergoes more and more 
oscillations (gyrations) before being reflected back. 

Clearly, for yp « t t/2 we can use adiabatic invariant 



theory (cf. Sec. |III C[ ) to treat the process of reflection in 
the mirror field in the throat of the cusp field, whereas for 
yo ~ the particle does not complete even one gyration 
about the magnetic field, so the adiabatic invariant is not 



defined on any part of the orbit. In order to determine 
on which part of any given orbit /i is approximately con- 
served, we compute the cyclotron frequency f c = w c /27r 
at each point on the orbit, where u) c (t) = \q\B(t)/m, and 
compare it with the time-rate-of-change of \nB. 

Suppose that the inequality f c > B/B holds over the 
interval t% < t < t% and is violated outside the interval. 
Adiabatic invariance theory applies during this interval, 
provided the particle has enough time to execute at least 
one gyroorbit. To determine the latter point, we calcu- 
late the total gyrophase change over the interval in which 
adiabatic invariance potentially applies, 



A0: 



■ dt 



(14) 



Then we define an adiabatic reflection as one for which 
A0/27T > 1 and a nonadiabatic reflection as one for which 
A^/2tt < 1. 





0.25 0.5 0.75 1 1.25 1.5 1.75 2 
FIG. 6. Adiabaticity parameter /S.<j>/2n vs. impact param- 
eter times 2/-7T for normally incident orbits with energy 0.01 
(solid line) and very low energies (dashed line). 

Figure ^ shows this adiabaticity parameter for the case 
of normal incidence, as depicted in Fig. ^. Two values 
of energy are shown, a relatively high energy £ — 0.01 
and the low-energy limit £ — ► (see below). Reflection is 
nonadiabatic for roughly 60% of particles in both cases. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 % 
FIG. 7. Adiabaticity parameter A<f>/2ir vs. impact param- 
eter times 2 /it for orbits incident at 20° to the normal with 
energy 0.01 (solid line) and very low energies (dashed line). 

Figure shows the dependence of A0/27T on impact 
parameter yo for particles incident at an angle of 20° to 
the normal in the x-y plane, and with p z = 0. (Here yo 
is defined such that the orbits have initial values x = Xo, 
y = Yq + yo, where X is an arbitrarily large negative 



6 



constant and the constant Yq is chosen so that yo = 
corresponds to an orbit symmetric about the x-axis.) 

It is seen that the adiabatic region is much reduced 
in the low-energy case, and has virtually disappeared in 
the high-energy case. At angles of incidence greater than 
25°, both high- and low-energy particles reflect nonadia- 
batically for all impact parameters. The ratio of the solid 
angle occupied by the cone of angles of incidence 77 < 25° 
to the cone of all possible angles of incidence, n < 90° 
is about 0.093. Thus we conclude that considerably less 
than 10% of particles reflect adiabatically. 

The low-energy limit referred to above is defined by 
observing that the lower the incident energy, the larger 
the value of — Re£ at which the particle reflects. Thus, 
in this limit we can replace sech (n£/a) by 2 exp(?i£/a) in 
Eq. (|^) and define the low-energy approximation as the 
result of replacing '5 with 

* low = 2V> x exp(^ , (15) 

where is defined in Eq. (^). 

The dynamics in this limit exhibits an important scale 
invariance: if we displace the orbit in the cc-direction us- 
ing the transformation £ = £' + h, where h is a real con- 
stant, then the flux function changes by a constant fac- 
tor: *iow(£) = exp(n/i/a)* low (£'). Inspecting Eq. |[(]) 
we see that the transformation p = exp(nh/a)p', t = 
cxp(—nh/a)t' leaves the form of the equations of mo- 
tion invariant. The energy is transformed according to 
£ — exp(2nh/a)£' . Clearly yo is invariant and it is also 
easy to show from Eq. ( )l4| ) that Atfi is invariant under 
this scaling transformation. Thus we have the result that 
in the low-energy limit the function A(/)(yo) is indepen- 
dent of energy. 

2. z-Motion 

Although our idealized system is infinitely long, any 
real system will be of finite length and it is therefore of 
interest to enquire as to the rate of drift in the z-direction. 
Figure ^ shows the motion in z for the case shown above. 
We see that, for the case p z = 0, the motion appears to 
be a random walk with no secular drift. 

z 



0.5 




FIG. 8. The z-component of the orbit shown in Fig. 



3. Periodic orbits 

As for the four-cusp system of Dahlqvist and Russberg 
p7| , we show that our n = 6 example is not completely 
chaotic by showing numerically that there is at least one 
stable orbit surrounded by invariant tori (Kolmogorov- 
Arnol'd-Moser or KAM surfaces) which make a small 
region around the periodic orbit dynamically inaccessible 
to the chaotic orbit filling most of the energy surface. 

We expect the most stable orbit to be the one with 
the highest symmetry allowed by the system, i.e. 2n-fold 
symmetry, since this is the smoothest orbit, least like the 
trajectory of a billiard ball. This is illustrated in Fig. || 
for the case n — 6, p z — 0, with the orbit passing through 
{x,y)o = (3,0). The corresponding momentum required 
to close the orbit is (p x ,Py)o = (0,0.04925), giving an 
energy £ — 0.00170082 and a period 28.96. Because of 
its 12- fold symmetry we call this the dodecagonal orbit. 



y 




FIG. 9. Stable periodic orbit with 12-fold symmetry — the 
"dodecagonal" orbit — as described in the text. The thick 
line shows the boundary of the energetically accessible region. 

To investigate the stability of such an orbit we linearize 
about the orbit and calculate the eigenvalues of the ma- 
trix evolving a neighborhood of the phase-space point 
(xo,p y o) in the x > 0, y = half plane to its intersection 
with the next surface of section that is equivalent under 
the symmetry operation £ 1— > exp(i7r/6)C- [The crossing 
time is found numerically by searching for the first zero 
of argexp(— i7r/6)£(t).] For the case shown in Fig. 0, the 
eigenvalues are —0.54162 ± 0.840624i, which lie on the 
unit circle in the complex plane, indicating stability. 

This stability is confirmed more graphically by the 
Poincare plots in Fig. |l0| for some neighboring orbits 
on the same energy surface as the dodecagonal orbit 
shown in Fig. [| Figure [l(] shows iterates of the map 
(x,p x ) 1— ► (x',p' x ) found by calculating the crossing time 
t with the 9 = tt/6 line as described above, then calculat- 
ing x' — Re exp(— i7r/6)C(i)i p'x — ^ e ex P( — i n /6)Pt(t)- 

Figure [l0| shows an island of regular motion in a vast 
chaotic sea — if we start much beyond the last orbit 
shown, the orbit rapidly moves far away from the peri- 
odic orbit. For example, the chaotic orbit shown in the 
Poincare plot, Fig. |], started at {x,y)o = (3.02,0), with 
p x = and p y adjusted to give the same energy as that 
of the periodic orbit plotted in Fig. [|. 
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The orbits in Fig. [L0| appear to lie on invariant curves 
that are topologically circular, being the intersection of 
invariant tori with the surface of section. The quasi- 
triangular shape of the outer orbits is due to the exis- 
tence of three unstable periodic X-points which define 
the separatrix between regular and chaotic motion (cf. 
the bifurcation with branching ratio 1/3 in Fig. 1(c) of 

H)- 
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FIG. 10. Poincare section transverse to the dodecagonal 
orbit showing several satellite orbits which appear to be on 
invariant tori. 

Scanning through a range of initial values of x, so as 
to vary the energy £, we find that for all £ less than 
the escape energy £ csc (and somewhat beyond) the do- 
decagonal orbit is stable. We have not done an exhaustive 
study of other periodic orbits, but have established that 
the 4-fold-symmetric "square" orbit is unstable below an 
energy of about 0.1, but stabilizes above this value. 

The existence of a stable periodic orbit shows that 
the multiple-cusp confinement system is not completely 
chaotic. However, the islands around the few stable or- 
bits are very small. For instance the area occupied in 
the x — p x plane by the island shown in Fig. ^ is about 
7 x 10 -5 , compared with the energetically accessible area 
of the Poincare section, J[2m£ — q 2 ip 2 (x, 0)] 1 / 2 ) dx « 
0.37, i.e. four orders of magnitude smaller. Even a very 
small amount of extrinsic stochasticity (e.g. small-angle 
collisions with other particles) will easily destroy such 
small islands of stability. 

Thus we conclude that, for practical purposes, the as- 
sumption of complete chaos in the free-particle dynamics 
is well justified, and hence assume that any orbit covers 
its energy surface ergodically. 

IV. WAVE-PARTICLE INTERACTION MODEL 

For a given wavelength, A, of the incident wave, the 
ratio \/{a/n) tends to infinity as n — > oo. Hence we 
consider the long- wavelength limit, in which the wave- 
particle interaction is via the uniform, oscillatory electric 
field E = Re (Eexpiwi), where E is a constant complex 
vector. On the other hand we assume u> 3> w P)e , where 
uj PtC is the electron plasma frequency, so that the oscil- 
latory part of the electrostatic potential can be ignored. 
Thus the electric field is taken to be produced by the 
vector potential Ai = Re [i(E/w) expiuit]. 

In reality, uniformity of E applies only during one wall- 
scattering event, which occurs over the scale length a/n, 



whereas E can be different at different points on the 
picket fence if A is comparable with a. In this paper we 
consider a model system in which E is strictly constant 
in space, but we model the real situation by choosing E 
from a random ensemble of values with probability distri- 
bution matching the distribution of field amplitudes and 
phases actually encountered on the cylinder r — a. We 
also assume the distribution of initial phases to be uni- 
form, which means that the ensemble averaging operator 
(■) automatically includes phase averaging. 

The model Hamiltonian determining the full particle 
dynamics is thus taken to be Hq + Hi , where Hq is given 
by Eq. (@), and 

Hi = — Re [i(p — qipe z ) • Eexpiwi] , (16) 

muj 

with e 2 being the unit vector in the z-direction. The 
above form for H follows simply by expanding (p — 
qA) 2 /2m and dropping the term quadratic in E as it 
has no spatial dependence and thus does not affect the 
dynamics. (If we had not taken E to be constant in space, 
it would be necessary to retain this term to include pon- 
deromotive force effects.) 

Note that H\ produces an oscillatory correction, 
dHi/dp, to the velocity, x, that is independent of po- 
sition, whereas the oscillation in p is localized around 
the edge region due to the localization of Thus p 
is nonoscillatory in the middle region where ip is negli- 
gible. This is a consequence of our choice of gauge for 
representing E, which automatically gives an oscillation- 
center representation [p3| , p. 47] in generalized momen- 
tum space. The localization of p to the edge region where 
the particles are reflected is desirable since this is the re- 
gion where irreversibility is introduced — in the middle 
region the particles simply respond adiabatically to the 
high-frequency field. 

It is possible also to remove the oscillation in the gen- 
eralized position coordinate in the middle region, so as to 
completely localize the interaction to the edge region, by 
making a canonical transformation to oscillation-center 
variables (^6). However, as we are interested in the dif- 
fusion in energy, not position, we have not found this 
transformation to be useful. 



V. QUASILINEAR DIFFUSION 

The Vlasov equation for the single-particle distribution 
function /(x, p, t) can be written 

d t f + {f,H} = 0, (17) 

where the Poisson bracket {/, H} = d x f-d p H—d p f-d x H. 

We decompose / into a nonfluctuating part, /, and a 
fluctuating remainder /. We define / = ((/)) where ((•)) 
includes not only an average over the wave phase via en- 
semble averaging, but a phase-space average over the un- 
perturbed energy surface: for arbitrary phase-space func- 
tion g(x, p, t), varying on both the fast timescale w _1 and 
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a slower (diffusion) timescale, we define ((g))(£,p z , t), 
varying only on the slow timescale, by 

((g)) = ^J^d 2 x jd 2 p5{£ -ff (x,p))( ff (x,p,t)) , (18) 

where r = (x 2 + y 2 ) 1 ^ 2 , d 2 x = dxdy, d 2 p = dp x dp y , 8 
denotes the Dirac ^-function, and the normalizing factor 
Af(£,p z ) is defined by 



AT = 



d 2 x ld z pS(£ - H (x,p)) 



(19) 



The result of applying this projection operation is a 
function only of £ and t, so that {/, Ho} = and hence 
the averaged part of the distribution function, /, is in- 
variant under the unperturbed dynamics. The projection 
of the distribution function onto the energy surface us- 
ing the averaging operator is an extreme form of phase- 
space coarse-graining. Owing to the highly chaotic na- 
ture of the unperturbed dynamics we assume that this 
coarse-grained distribution function relaxes to a function 
of the constants of the unperturbed motion, £ and p z , 
on a timescale much faster than the quasilinear diffu- 
sion timescale. We assume all particles to be confined, 
so that the region of phase space defined by £ = const, 
p z = const is bounded within r < a. 

Applying the operation ((•)) to Eq. (|l7|) we have 



dt/+«{/,ffi}» = 



(20) 



Writing the Poisson bracket in the form {/, Hi} = 
<9 X • (fdpHi) — d p ■ (fd x Hi) and integrating by parts (as- 
suming the particle confinement is good enough that the 
boundary contribution can be ignored) we find 



(({lH!})) = ~ (//((ft)) 



(21) 



where £ is the rate of change in the energy integral of 
the unperturbed system, £(t) = H n (x(t), p(i)), along the 
perturbed orbit. Noting that H = {H ,H + Hi} = 
{H ,Hi}, we see that 



£ = {H ,Hi} 



(22) 



Also, p z = {p z , Ho + H\} = \p Zl Hi \ [which vanishes for 
our simple interaction term, Eq. ( |16|) 1 . 

Subtracting Eq. (I2O) from Eq. (|l7|) we also have 



d t f + {f,H } = -{f,Hi} + O(V 2 



(23) 



Linearizing Eq. (|23|) and solving by integration along the 
unperturbed trajectories from an initial time — T, we 
have 



/( X ,p,t)=/(x(-T),p(-T),-T) 



d£ 



—T 



£'^(£', P 'z,t') 



+p z — (£, Pz ,t) 



(24) 



In calculating £' = £(t') using Eq. (P4), the right hand 
side is to be evaluated at the point (x(t'),p(t')) on the 
unperturbed phase-space trajectory that passes through 
(x, p) at time t, and similarly for p z (t') if an interaction 
model is used for which it is nonvanishing. 

We now observe that, for large T, J(x', p', t')(t' = -T) 
becomes decorrelated from £(t) and p z (t) and thus does 
not contribute to the averages on the right hand side of 
Eq. (21). The decorrelation time is the duration of one 
wall-scattering event, which is of the order of the transit 
time 



Ttr(£) = -r-r 
n\v\ 



(25) 



of a free particle with speed \v\ = (2£/m) 1 / 2 through the 
scale length a/n of the magnetic field variation. Thus, 
assuming T 3> r tr , we can set the initial value term 
f(x',p',t')(t' = —T) to zero without significant error. 

Also, if x(i) is in the wall-interaction region, where £ 
and p z are significant, then x(— T) is far from the wall so 
£{—T) and pz(— T) are negligible (because ip is essentially 
zero there — see the discussion at the end of Sec. IV). 
Thus we can, to a very good approximation, extend the 
lower limit of the integral in Eq. (|24|) to —00. 

Whereas T is assumed large with respect to Tt r , we 
assume it to be small with respect to the characteristic 
evolution time for the distribution function /. (That is, 
we assume the wave to be of sufficiently low amplitude 
that it takes many wall-interaction events for significant 
heating to occur.) Thus we can also make the Markovian 
approximation that f(£',p' z) t') can be moved outside the 
integral in Eq. ( ^4| ) with negligible error. 

Substituting Eq. (§|) in Eq. @ and then Eq. @ we 
find (assuming p z = 0) the quasilinear diffusion equation 



dt 



_L__d_ 

Jfd£ 



ND 



df 
d£ 



(26) 



where D(£,p z ) is the energy diffusion coefficient, defined 
by 



C(r)dT 



(27) 



with the two-time correlation function C(£,p z ,r) 

C(r) = ((£(t - r)£(t))) = ((£(r)£(0))) , (28) 

where the second form follows from the fact that, because 
of the average stationarity of the dynamical system, C 
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depends only on the time difference, r = t — t' . Also 
note that the projection operation ((■)), Eq. (|l8|), can 
be done using either initial or final values as indepen- 
dent variables because the Jacobian of the transforma- 
tion is unity (preservation of phase space volume) and 
Hq is an invariant of the unperturbed motion. Thus 
C(t) is an even function, which fact we used to ex- 
tend the upper limit of the integral in Eq. ( p7| ) to in- 
finity. We can also use time reversal invariance to show 
C(£,p z ,r) = C(£, - Pz , -t) = C(£, - Pz ,t). 

We end this section by calculating the heating rate due 
to Fermi acceleration. First we define the total plasma 
energy per unit length in the z-direction, 



U(t) 



d£ 



d Pz Af£f . 



(29) 



Differentiating U with respect to time, using Eq. (|26|) and 
integration by parts we find the rate of power deposition 
into the plasma due to reflections from the confining edge 
magnetic field under the influence of an electromagnetic 
wave 



U{t) 



d£ / dp z ND 



0£ 

d£ 



(30) 



We have assumed that df jd£ vanishes at an energy less 
than or equal to t he m aximum confined energy q 2 i\}\j2ra 
discussed in Sec. II C so that we can ignore boundary 
contributions. 



VI. ONE-DIMENSIONAL MODEL 



We saw in Sec. [II D 1 that most particles reflect nona- 



diabatically in less than one gyroperiod, and thus should 
not be sensitive to the details of the y-variation of ip 
(i.e. subtle resonance effects should not be important for 
most particles). This suggests we estimate the effect of 
nonadiabatic reflection by using a one-dimensional model 
Hamiltonian obtained by replacing ip(r, 8) in Eq. (|l|) 
with an axisymmetric flux function, ip(r) (cf. the one- 
dimensional model used by Yoshida et al. Q ) . The con- 
servation of the angular momentum pg then allows formal 
integration of the equations of motion by the method of 
quadratures. 

Although such a one-dimensional flux function violates 
Laplace's equation, and therefore would require a plasma 
current to produce it, this fact is irrelevant to the single- 
particle dynamics. By suitable choice of ip(r) we can 
model the gross radial confinement properties of the two- 
dimensional flux function. The main loss in the physics 
is that ijj{r) produces no radial component of B, and 
hence no interaction with Eg . But if we assume perfectly 
conducting wall boundary conditions, E must be purely 
radial at the vacuum vessel wall (assumed just inside the 
array of magnetic dipoles) so Eg (and E z ) can be assumed 
to be small in the interaction region anyway. 



We further simplify the reflection dynamics by going 
to the large- n limit, so that the dipoles become a lin- 
ear array a nd we can use Cartesian coordinates as in 
Sec. [Ill D l[ Also, Figs. | and indicate that the low- 
energy approximation, Eq. ([l5|), is good for nonadiabatic 
reflections. In this limit the field strength B is indepen- 
dent of y, so we define the equivalent one-dimensional flux 
function ip(x) as that which gives the same field strength 
B(x). That is, B = |*( ow (0l = &{x), where 



/ TLX \ 

ip(x) = 2?/>x exp^— J 



(31) 



(In the above we have shifted the origin of the x-axis to 
lie on the same line as the array of dipoles.) As a final 
simplification of the unperturbed dynamics we evaluate 
D only at p z = 0. That is, we consider only unperturbed 
orbits having the constant of the motion p z = 0. 

In the large-n limit the boundary region where ip is not 
small makes negligible contribution to Eq. (|l^) and thus 
we find the normalizing factor M to be independent of 
energy, 



M=2n 2 



(32) 



The equations of motion Eq. (|1C|) can be integrated 
explicitly to give 



x(t) = (— \ In ( u sech U ^ — ^ m 
V n J V tx 



(33) 



Px{t) = _ 2w ( ™) tanh 2 ^ - hue! (34) 

\nTxJ t x 

where rx is defined in Eq. (^) and the constants of in- 
tegration are u = exp(na; max /a) and i max - Inspecting 



Eq. 



we see that x r 



is the maximum value of x 



attained over the entire orbit, and i max is the time at 
which this point is reached. 

Assuming a perfectly conducting vacuum vessel we 
set E y = E z = 0. Then, from Eq. @, H 1 = 
— (q/muj)Ke(ip x E x expiu!t) and we have the simple ex- 
pression for the instantaneous power transfer to a parti- 
cle, Eq. P), 



a . dHi q 
£ = -Px^ — = — Rc 

Op x LU 



iE x x(t) exp iwt 



(35) 



In evaluating the diffusion coefficient using Eq. ( |27| ) , it 
is convenient first to commute the time integration with 
the averaging operation, so that we first consider the time 
integral of £, which gives the total energy change, A£ , 
in one collision with the wall. Inserting the analytical 
solution Eq. (33) in Eq. ( |35| ) and integrating from t = 
—go to +oo, we get 



A£ 



/7TWTX\ 

cosech 

\ 4u / 



x Re 



(iqE x 



exp iuit Ti 



(36) 



10 



Since Eq. (Q) expresses the orbit in terms of the con- 
stants of integration rather than the initial conditions, 
to evaluate the phase-space average we change variables 
from the initial conditions x,p x to u and s = uimax/Tx so 
x = (a/n) ln(usech 2s), p x — 2u(ma/nTx) tanh 2s. The 
Jacobian of this transformation is 4(ma 2 /n 2 Tx), so, us- 
ing Eq. ( |32"| ) in Eq. (18), the phase space average over 
wall scattering events is transformed to 



((•» 



8 / ma 



it \n*Tx 



duO(£- Au 2 £ csc ) 
/2m~(£ -4u 2 £ esc ) 



1/2 



ds 



(37) 



where £ esc is defined in Eq. (0) and 6(-) is the Heaviside 
step function. 

Using Eqs. @ and (§71) in Eq. @ we have 



D = 



3tt 



TTCJTtr 



(38) 



where \v\ = (2£/m) 1 / 2 is the mean velocity in the field- 
free region and r tr (£) is the transit time defined in 
Eq. ( f25| ) . (Note that D does not depend on the strength 
of the magnetic field in this model, only the scale length, 
because a change of ipx is equivalent simply to a shift in 
the origin of the x-axis by an amount of order a/n.) 




0.5 1 1.5 2 2.5 3 
FIG. 11. Function G(w) denned in the text (solid line) and 
the large-ui asymptotic approximation (dashed line). 

The function G is defined as a one-dimensional inte- 
gral, 



G(w) 



3w 2 
2 



1 p cosech 2 (w I p) 



(1 



p 2y/2 



dp 



(39) 



and is plotted as the solid line in Fig. O. The func- 
tion has been defined so as to approach unity as w — > 0, 
as discussed below in the context of the low-frequency 
limit, lo <C 1/Ttr- The asymptotic behavior shown by 
the dashed line is discussed below in the context of the 
high-frequency limit lo 3> 1/T tr . 



A. Low-frequency (Fermi) limit 

Fermi was concerned with the collision of cosmic 
rays with relatively slowly moving gas clouds. In our 
problem this corresponds to the low- frequency limit, in 



which a particle scatters off the magnetic field in a time 
much less than the period of the applied field. This makes 
the argument of G, w = nujTtr/2, small. 

For \w\ <C 1, we can approximate cosech 2 (w/p) in 
Eq. (B9) by p 2 /w 2 over nearly the full interval. Eval- 
uating the integral we find G(0) = 1. 

Thus, in the low- frequency limit, 



D = (4/3ir)(q 2 (\E x \ 2 )/u J 2 )(\v\ 3 /a) 



(40) 



This result may be understood as the Fokker-Planck dif- 
fusion coefficient ((A£) 2 )/2t co h for particles oscillating 
in the applied field with the quiver velocity v, given by 



i(\e x ?) 1/2 



(41) 



giving the typical energy step at each collision with the 
wall AS = m\v\v. Taking as a typical time between wall 
collisions r co n = a/\v\ we recover, to within a factor of 
order unity, the low-frequency energy diffusion coefficient 
above. 



B. High-frequency limit 

In the high-frequency limit, the particle oscillates many 
times during a collision with the magnetic field and we 
would expect it to respond to the applied field adiabati- 
cally, gaining little energy. 

For \w\ 3> 1 the dominant contribution to the in- 
tegral comes from a narrow boundary layer near p — 
1, in which cosech 2 (w/p) may be approximated by 
cxp(— 2\w\) exp(— 2\w\(p— 1)). This gives the asymptotic 
behavior 



G(w) 



3V^H 3/2 exp( 



-2H) 



(42) 



From the dashed line in Fig. [n]we see that this asymp- 
totic formula gives good agreement with the numerically 
calculated result for \w\ greater than about 1. We see 
that the energy diffusion is indeed exponentially small in 
this limit. 



VII. DISCUSSION 

In this section we give the magnetic parameters of the 
theory for a typical experimental device and make some 
observations as to the possible implications of the theory 
for such experiments. 

Multipolar magnetic cusp confinement has become a 
conventional method for reducing plasma loss on the 
chamber walls and keeping the inner plasma volume free 
from magnetic field |13|. This was used in the electron- 
cyclotron resonance (ECR) plasma formation experiment 
ECRIN (ECR Ions Negatifs) In ECRIN, mi- 

crowave argon and hydrogen plasmas were excited in a 
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cylindrical vessel of internal radius of about 6 cm and 
length 17 cm surrounded by twelve radially magnetized 
linear permanent magnets of alternating polarity, max- 
imum magnetic field strength 0.2 T and microwave cw 
power of 100-1000 W at a frequency of 2.45 Ghz was 
delivered at one end of the vessel. 

The primary heating occurred near the microwave in- 
put window, but it is of interest to consider whether col- 
lisionless secondary heating of free particles is possible 
further down the tube, which we can model by idealiz- 
ing the permanent magnets as the n — 6 linear magnetic 
dipole configuration used for illustration in the present 
paper in Figs. |l|-^ and Figs. p}-|lO|. Using a = 6 cm 



gives the length unit in these figures (see Sec. 1IC) as 
a/n = 1 cm. (In this paper we have ignored collective ef- 
fects, collisions and atomic processes, all of which may be 
important in such experiments, so the use of the ECRIN 
parameters should be regarded as illustrative only.) 

For ECRIN, the magnetic dipole strength was esti- 
mated to be around K = 1.5 x 10~ 5 Tm 2 . Using this 
value in Eq. (Q) gives the escape energy for electrons as 
^csc ~ 198 keV, and that for singly charged argon ions as 
£L r « 2.7 eV. 



For electrons of energy 5 eV the transit time, Eq. (25), 
is r-tr ~ 7.5 ns, so that for the microwave heating fre- 
quency of 2.45 Ghz the argument of the transit-time re- 
duction factor G in Eq. (pq ) is w = nujTt r /2 i=s 182. This 
gives G(w) w 5.5 x 10~ ! Thus we see that nonreso- 
nant Fermi acceleration is clearly not an important ef- 
fect in such ECR experiments. On the other hand, with 
l/r tr « 133 MHz, this effect can be important in rf heat- 
ing experiments. 

Given the strong transit-time suppression of nonreso- 
nant heating, it may be of interest to consider the reso- 
nant heating of the few free electrons penetrating deeply 
enough into the cusps to reach the ECR layer, and this 
could in principle be calculated using the quasilinear for- 
malism developed in this paper. 

However we shall content ourselves here simply with 
estimating the proportion of the ECR layer that is ac- 
cessible to the free electron s, as opposed to the trapped 
electrons discussed in Sec. Ill C . In the neighborhood 
of the ECR region, Eq. (|l^) is satisfied only in the nar- 
row cusp regions directed toward the magnets. We can 
thus Taylor expand %j) to approximate this inequality by 
r\A8\uj c (r) < 2(2m£ /to) 1 / 2 in polar coordinates, where 
uj c = \q\B/m is the electron cyclotron frequency (= to in 
the ECR layer). 

Summing these angular ranges over all the 2n cusps 
and dividing by 2tt gives the fraction of the ECR layer 
accessible to free particles. Approximating r by a gives 
this fraction to be (4/7r)(wTt r ) _1 « 1%. On this basis we 
would expect nearly all the ECR power to be deposited in 
the trapped particles, with the free particles being heated 
through heat conduction from the trapped population 
and other such indirect processes. 

This paper has focussed only on the effect of chaos 
as the source of stochasticization. In reality, particle- 



particle collisions may be equally or more important. Our 
collisionless energy diffusion coefficient will still be valid 
as an additive contribution to the total energy diffusion 
coefficient provided A m f p S> a/n, for then most parti- 
cles transit the high-field edge region without suffering a 
collision. Elastic collisions in the central region simply 
provide a further stochasticization and do not affect our 
result provided the above inequality is satisfied. Rare 
collisions within the edge region would provide an inde- 
pendent additive mechanism for energy diffusion which 
might or might not dominate our collisionless mechanism 
depending on the ratio of transit time to the period of 
the applied wave. 



VIII. CONCLUSION 

We have shown that, in such strongly nonaxisymmet- 
ric experiments as the multicusp geometry analyzed here, 
there is a strong collisionless stochasticization process 
due to the chaotic nature of the unperturbed particle mo- 
tion. This justifies the use of the random phase approx- 
imation for successive kicks produced by coherent wave- 
particle interaction without having to invoke a nonlinear 
threshold for resonance overlap, or collisions. Such sys- 
tems cannot be analyzed by area-preserving maps, and 
thus fall outside the general framework usually assumed 
for the analysis of rf and microwave heating in bounded 
systems p2[ . 

As an alternative to the Fokker-Planck approach for 
deriving the energy diffusion equation we have developed 
a variant of the quasilinear diffusion formalism based 
on averaging the single-particle Liouville equation. This 
provides a general and efficient formalism for treating 
complex geometries. 

We have applied the formalism to an exactly soluble 
model for nonresonant Fermi acceleration and found a 
transit-time correction factor that becomes exponentially 
small in the high-frequency limit. 

Finally, we have illustrated these concepts using pa- 
rameters from a fairly typical electron-cyclotron heating 
experiment. 
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